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A LIMIT THEOREM FOR THE SUM OF SQUARED DIFFER- 
ENCES OF AN INTEGRATED ITO PROCESS WITH APPLICA- 
TION TO INVERSE SCATTERING 

JOHN F. A. FLETCHER,* University of Cambridge 

I 

Abstract 

We investigate a functional obtained by summing the squared differences of 
the integral of an Ito process over disjoint intervals. The limit of this sum is 
shown to converge in probability to two thirds the quadratic variation of the 
underlying process. An application to inverse scattering from a random fractal 
surface is presented. 
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Introduction 

Imagine that the real-valued function (A" t ) t6 R represents the density of a physical 
quantity. For example, X could be the density of a fluid along a one-dimensional 
subspace parametrised by s, or the energy density of radiation in a scattering problem. 
With regard to measurements, such densities are, of course, idealisations. In reality, 
one would attempt to estimate the value of X, for fixed s, by taking some sort of 
local average. Define a discrete set of times % = Afc where k = 0, 1, 2, and A > 0. 
Denote the average of X over the k th interval \tk,tk+i) by Xk where 

X k :=^ J X s ds. (1) 
tk 
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1. Local average and quadratic variation 

Quadratic variation is a useful notion in many circumstances, for example, it can 
provide a way to estimate the volatility of the quantity under study. If local averages 
of the type ([T]) are available, the following theorem may be applied: 

Theorem 1.1. For a standard Ito process (A^) f >o where 

X t =X Q + f a s ds+ f b s AB S (2) 



o Jo 



with average X k defined by (flj we have 
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Y^iXk-Xk^) 2 ^ -(X) t as A^O (3) 



fc=i 



where {t k ;k = 0,1,..., n = t/A} specifies a partition of [0,t], and — > denotes conver- 
gence in probability. 



Proof. Integrating X by parts: 

tk+l tk+1 



J Xsds^sX^T 1 - J sdX s 



(4) 



It follows that 



tk tk + l 

X k -X k _ 1 = ^ J (s-t k -i)dX s + ^ J (t k+1 -s)dX s 

tk-l tk 



tk+l 

Ps (k,A) dX s 



tk- 



whcrc 



(s-t fe _i)/A if t k -i < s < t k , 
(t k+i -s)/A \it k < s <t k+1 . 



p s (k,A) = 

Assume for definiteness that n is even (the case of odd n follows in an obvious way) . 



Squared differences of an integrated Ito process 



Splitting the sum in ([3]) into two: 



n-l 



53 ("^ fc — Xk-l) 2 
1 

n-l ( tk + 1 

53 J Ps(k,A)dX s 



k=l 
n-l 



k=l . t 
n/2 I * 



n/2-1 / t2fc +! 



53 J p s (2k-l,A)dX s \ + \ J Ps-A(2k-1,A) 

k—1 \ + _, _ / k— 1 \ +„, 



dX., 



(5) 



Define the continuous function r] t (A) as 

oc 

Vt (A) = l {t2k _ 2 <t<t 2k }Pt(2k - 1, A) (6) 

k=l 

and the process X t := X + J * Vs(A) dX s . Also, we use the following notation 



Qa (Xt) '■— 51 \ X kA - -X"(fe-1)A 



A' 



(7) 



fc=i 



Set Ufe := ^2fe- The first sum in ([5]) becomes 

n/2 / \ 

X] / p s (2fc-l,A)dX s = ^ 
fe=i \ , ' I 

\ t 2k — 2 J 



n/2 



k=l 



(8) 



The the limit A — >• of (J5J looks a lot like a quadratic variation, except that X is 
dependent on the grid spacing A. 

The process X can be expressed as X t = A t + Z t where A has finite variation and 
Z is a local martingale. Clearly \Z Uk — Z Uk _ 1 \ < \Z Uk — Z Ukl \ and \A Uk — A Uk l \ < 
\A Uk — A Uk _ 1 1 for all k. As one might expect, given the well known properties of the 
quadratic variation, the process A does not contribute to ((8|) in the limit A — > 0. This 
can be seen in the following waj: 



?l(Zuk Z Uk _ 1 )(A Uk — A»»_i, 



< [sap\Z Uk -Z Uk _A )-Q$i(A t ) 



< I SUP \Z V ^ - Z; 



".. ill' Q2A (^*) 



Compare proposition 1.18, p. 128, [3]. 
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and this converges to zero as A — > because of the continuity of Z . Similarly 

hm^(l Ufc ) 2 = 0. (9) 

k 

We now require two well known results from the literature, which are stated as the 
following Lemmas: 

Lemma 1.1. Let (M s )o< s <t be a uniformly continuous martingale such that \M S \ < 
K < oo, P-a.s. then 

hmE[Qi 4) (M t )]=0. 
Proof. See for example Lemma 5.10, p. 33, [1]. 
Lemma 1.2. Let (M) s < K for all s E [0,t], P-a.s. then 



E 



(Qi 2) (M t )-(M> 



2 



< 2E[Q { l\M t )] + 2E[(M) t ■ m t (A, (M)) ' 



where m u (A,X) :— su~p{\X t — X s \;0 < s < t < u and \t — s\ < A} is the modulus of 
continuity. 

Proof. See proof of theorem 5.8, p. 34, pQ. 
Assume that Z is a bounded, uniformly continuous martingale. Lemma 1.1 implies 

E[Q { *l(Z t )] < E[Q$(Z t )] ^0 as A -> 0. (10) 

Furthermore 

m t (2A, (Z)) < m t (2A, (Z)) — > as A^O (11) 

by the sample path uniform continuity of (Z). Assume that {Z) is bounded, then so 
is (Z), and, since TO((2A, (Z)) is bounded 

limE\(Z)fm t (2A,(Z})] =0 (12) 

A— >0 

by the dominated convergence theorem and ([TTj) . Putting these facts together and 
using Lemma 1.2 we have 

2" 



lim E 

A->0 



(Qfl(Z t ) (Z) 
Convergence in L 2 implies convergence in probability. 



0. (13) 
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Next, we would like to show that (Z) converges, in some sense, to a well defined 
limit, independent of A. 

n/2 y 



n/2 u k n/2 u p 



£ / ry s 2 (A)(^-^)d S + ^ J V 2 s (^b 2 K ds (14) 



fe=i „ fc=i „ 



where u£ £ [wfc-i, U&) is such that b u * = sup{fo s ; Uk~i < s < Wfe}. The first integral on 
the right hand side of (|14p has the bound 

n/2 «* n/2 u * 



1 - ' J ' k=i „ 

/ & 2 d S -^2A&i 
Jo k=i 
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2 a ' S; as A^O. 



For the second integral in (|T3)> we have 

n/2 u * 



2 / r, 2 (A)6 2 ,d S (15) 

i 11/2 1 /■* 

- V 2A& 2 , ^i- b 2 s ds as A -> 
3^ B * 3i s 



and so 



g^(Z t )^-(^) t as A^O. (16) 

The technique of localisation can be used to generalise (TTBl to the case where Z and 
(Z) are not necessarily bounded, using the stopping time 

Tj = M{t > 0; \Z t \ > j or (Z) t > j}. (17) 

The limit of the second sum in (0 can be evaluated in precisely the same way as the 
first, and converges to the same limit. Finally © follows. 

2. Application to inverse scattering 

Norris [2J showed that, in the ray-theory limit, fluctuations in the intensity of 
light scattered from a corrugated surface with Brownian motion slope (crW s ) s eWL are 
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described by the Ito process (i?t)teR where 



dRt = — (1 - RAdt + -jRtdB t 




(18) 



and a is a positive constant. Measurements of the average energy flux between points 
{tk',u = to,ti, ...,t n = v}, such that tk — ffe-i = A > 0, correspond to Rf.. The 
volatility of the surface slope can be obtained approximately for small A, from (J3) 
since 
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